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ABSTRACT 

We discuss the possible origin of the radial tr uncations in stellar discs, using measure - 
ments that we presented in an earlier paper ijKregel. van der Kruit fc de Griisl l2002). 
A tentative correlation is found with the de-projected face-on central surface bright- 
ness; lower surface brightness discs tend to have a smaller truncation radius in units 
of scalelength. This and our earlier finding that in smaller spirals the truncation tends 
to occur at a larger number of scalelengths are best reproduced when the truncation 
is caused by a constant gas density threshhold on star formation. 
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1 INTRODUCTION 

Stellar discs have a finite size. In the extreme outer 
parts the stellar light distribution diminishes more steeply 
than the exponential decline over the main disc and 
drops to low values beyond the s o-called truncation ra- 
dius -Rma. v llvan der Kruid Il978l . |vaii^ej^<rui^&Searld 
1981allbl Il982t hereaf t er KS1-3; iPohlen et al.1 l2000t 
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20011. hereafter GKW: iFlorido et al.ll200ll) . This truncation 
or cut-off of the disc light is often directly visible in contour 
maps of edge-on spir als and usually occur s at a radius of 3 
to 5 disc scalelengths llvan der KrufifeOOlal) , although stellar 
discs that extend to a much lar ger number of scalele ngths 
are certainly known to exist (e.g JWeiner et al"ll2001b|) . The 
truncation is most easily found in edge-on spirals because of 
the line-of-sight projection and the associated higher surface 
brightness. 

In less inclined spirals, for which azimuthally aver- 
aged radial light profiles are routinely studied, the non- 
axisymmetric component (e.g. spiral structure, lopsided- 
ness) can smooth out a truncatio n present in t he old disc. 
This effect was first noted by Ivan der Kruitl 1^.988). For 
16 face-on spirals, of which 15 did not show any sign of 
a truncation in the azimuthally averaged light profile, he 
found that the three outermost isophotes were much more 
closely spaced than the inner ones, providing clear evi- 
dence for the truncation. The Milky Way disc probably 
also shows a truncation, with recent estimates of 7? ma x 
rangi ng from 10-15 kpc based on near-infrared star counts 
l)Ruphv et al.l fl9961 and the near- and far-infrared sky sur- 
vey of the COBE/DIRB E instrument jFreudenreicblll99St 
iDrimmel fc SpergellSool . Recently, wide-field surveys have 
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also r evealed trunca tions in the Loca l Group galaxies NG C 
2403 toavidedtooljl) and NGC 3109 iPemers et alJl2003ft . 



The origin of the truncation of the stellar disc is still 
unclear and c ould possibly result from a number of physi- 
cal scenarios llvan der KruidlioOlal : lFreemanll2002l) . For ex- 
ample, if there has been no major re-distribution of the 
disc angular momentum during its formation and evolution, 
the truncation may reflect the maximum specific a ngular 
momentum of the protogalaxy llvan der Kruid Il987l) . This 
woul d imply that the H I extending beyond the stellar disc 
(e.g. iBroeils fc Rhee!ll997D has been accreted, and that the 
truncation i s associated wi th a small drop in the rotation 
curve (KS3: lBottemalll996l) . Another interesting possibility 
is that the truncation corresponds to the radius at which 
the gas density dr ops below a threshold d ensit y necessary 
for star formation (iFall fc Efstathioulll980l: KS3; iKennicutd 
119891: ISchavell2004D . 



In iKreeel. van der Kruit fc de Griisl ll2002l . hereafter 
KKG) we fitted truncation radii to the photometry of 3 4 
edge-on galaxies from the sample of Ide Griisl <ll997t Il998l) . 
The selection criteria and properties of this sample are sum- 
marized in KKG (sect. 2 and table 1). The sample covers 
a large range in rotation velocity (w max = 50-400 km s -1 ) 
and is dominated by spirals of intermediate- to late-type. 
Of the 34 galaxies we found that at least 20 have truncated 
discs. The truncation radii -R max have been listed in Table 5 
in KKG. Here we further compare the observed truncation 
radii to those predicted by the scenarios proposed for the ori- 
gin of the truncation, and specifically address the behaviour 
with disc surface brightness (see also lKregelll2003t) . 
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Figure 1. (a) - -R ma x/^R versus disc scalelength. The arrows indicate lower limits for galaxies for which no i? m ax could be determined. 
The cross-hatched region shows the predicti on of the coll apse model for disc galaxy formation (see text). The dashed lines show the 
predicted threshold radii for star formation JSchavc 2004J) for disc masses of, from bottom to top, 7.5 X 10 9 , 5 X 7.5 10 9 and 25 X 
7.5 10 9 Mq. (b) - i?max//iR versus face-on central surface brightness in the /-band. The arrows and the cross-hatched region is as in 
(a). The dashed line shows the constant star-formation threshold prediction. The shifted dotted line is obtained when allowing for an 
underestimation of the surface brightnesses by 0.7 mag (see text). 



2 THE TRUNCATION RADIUS VERSUS DISC 
PARAMETERS 

We show the distribution of -Rmax/foi versus scalelength 
in Fig. [Ql. For those galaxies for which no truncation was 
found, the lower limit is shown (arrows). Figure^ appears 
to show a subtle increase in i? max /ftR towards small scale- 
lengths; the average ratio for the spirals with /ir < 4 kpc 
is 4.4, more than one standard deviation higher than that 
of the entire sample. This increase of i? max //iR may be re- 
lated to its decrea se found at very large scalelengths by 
IPohlen et alJ (UoOO) . However, the reality of this feature is 
not entirely clear considering the modest sample size and the 
selection effect against galaxies of small physical sizes and 
low surface brightness. We are, for example, still missing 
small low surface brightness galaxies, which may have en- 
tirely different i? ma x/^R. While our view of the distribution 
of -Rmax/ftR is certainly obscured by this selection effect, ad- 
ditional information can be obtained from face-on samples. 
We have estimated a l ower limit to f? ma x for the spir als 
in the face-on sample of Ide Jong fc van der Kruitl jl994l) . a 
sample which is dominated by spirals of small scalelength. 
By taking the lowest contour in their i?-band contour plots 
as a lower limit to i? m ax, we find an average (i? m ax/^R) min = 
4.0±1.1 (la). This lower limit and the tentative increase of 
J?max//iR towards small scalelength spirals suggest that the 
ratio of truncation radius to disc scalelength in small scale- 
length nearby spirals is at least four. 

Figure ^> shows the dimensionless truncation radius 
versus the de-projected face-on centr al surface brigh tness 
(corrected for Galactic extinction as in Ide Griislll99a) . The 
ratio -Rmax//iR appears to be correlated with surface bright- 
ness in the sense that lower surface brightness discs tend to 
truncate at a smaller number of scalelengths. A Spearman 
rank correlation test appears conclusive, yielding a corre- 
lation coefficient of 0.61 (or a confidence level greater than 
99%). However, the errors are substantial in many cases and 



not entirely random, so that the test could have given a false 
positive result. Note that the face-on surface brightnesses 
are probably underestimated, by ~ 0.7 mag according to 
a comparison betw een the Tully-Fisher relations of edge-on 
and face-on spirals jKreeel. van der Kruit fc Freemanll2004l ; 
hereafter KKF). The large lower limits for ESO 340-G08 and 
ESO 555-G36 are puzzling. The discs of these spirals appear 
to be far more extended than the norm. 

Figures & b combined suggest that in small scale- 
length, high surface brightness spirals the truncation occurs 
at at least four scalelengths. This is important, because small 
spirals (fiB, ^ 4 kpc) are th e most numerous in the lo cal Uni- 
verse <|van der Kruitl Il987t Ide Jong fc Lacevl 1200(1 KKG). 
As another corollary, the range in face-on surface brightness 
among spirals is narrower at -R max than at the disc centre. If 
the discs in the current sample are exponential out to J? ma x 
then the distribution of /-band face-on surface brightnesses 
at J? m ax has an average of 25.3 mag arcsec -2 and a la scatter 
of 0.6 mag arcsec -2 . This dispersion is considerably smaller 
than the 1.1 mag arcsec -2 dispersion in the central surface 
brightness. This also implies that in the face-on view, the 
truncation is just as easily (or laboriously) detected in LSB 
as in HSB spirals. 



3 COMPARISON WITH HYPOTHESES FOR 
THE ORIGIN OF THE TRUNCATIONS 

We can compare the observed trends to the pre- 
dictions of the analytical collap se mo del of disc 

galaxy form ation llFall fc Efstathioul Il98(l iGunnl Il982l : 
iDalcanton. Snergel fc Summers! Il997l . hereafter DSS97). If 
protogalaxies are relatively sharp-edged, then the collapse 
theor y also predicts the outermost radius of the baryonic 
disc (van der Kruit Il987l) . This radius corresponds to the 
material with the highest specific angular momentum in 
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the protogalaxy. To quantify this prediction we calculated 
model surface density profiles using the method of DSS97 
for a range in spin parameter and total mass of the proto- 
galaxy. We assume a constant baryonic mass fraction F = 
0.10 and a constant efficiency for turning the baryons into 
stars over the age of the disc at (M/L)di sc = 2 (see KKF). 
The scalelengths were derived u sing a method similar to 
the 'marking-the-disc' method (iFreemanl 11970). and the 
outermost radii were obtained by taking the radius at which 
the density drops to zero. The result is indicated by the 
cross-hatched region in Figs.0 

In the collapse theory, both the outermost radius and 
the scalelength of the baryonic proto-disc increase with 
the mass and angular momentum of the protogalaxy such 
that their ratio remains approximately constant at 3-4. 
Taking a different or non-constant baryon mass fraction 
does not significantly change this result. The predicted 
fl max //iR is sli g htly s maller than the ratio of 4.5 predicted by 
Ivan der Kruid l|l987T ) based on a comparison of the angular 
momentum distribution of an exponential disc with that of a 
uniformly rotating, uniform sphere with A = 0.07. Although 
the prediction of the collapse model roughly coincides with 
the average observed -Rmax/^R, it can not explain the exis- 
tence of discs which extend to a relatively large or a small 
number of scalelengths. In particular, it does not predict an 
increase of i? max //iR toward small scalelengths (Fig. or 
a decrease towards low surface brightness discs (Fig. . 
Taking a different halo density profile and/or angular mo- 
mentum profile could change the value of -R ma x//iR, but not 
its constancy. Including additional prescriptions for mass ac- 
cretion, star formation and supernova feedback can lead to 
Rmnx/hn ratios that do change as a functio n of central sur- 
face brightness, but with the opposite sign l|van den Boschl 
l200lF) . Perhaps taking a range in angular momentum pro- 
files, as suggested by N-body simulations in the ACDM cos- 
mology tBullo ck et all2001aF) . or re-distributing the angular 
momentum during and/or after the collapse could solve the 
discrepancy. 

Alternatively, the stellar disc truncation may be caused 
by the inhibition of widesp read star formation be low a 
critical gas surf ace density jFall fc Efstathioul ll980: KS3; 
lKennicuttlfl989l) . if the corresponding critical radius is ap- 
proximately constant over time. This star formation thresh- 
old is sus pected to be re lated to the stability of the gas disc. 
Recently, ISchavd a2004f) made a prediction for the threshold 
radius based on simulations of the thermal and ionization 
structure of the gaseous discs ass embled in the galaxy for- 
mation model of iMo et al.l il998h . In these simulations the 
transition to the cold ISM phase is responsible for the onset 
of local gravitational instability that triggers star formation. 
This transition to the cold phase is independent of the shape 
of the rotation curve and occurs at a critical gas surface 
density, which for reasonable values of the gas fraction, tur- 
bulence, metallicity and the UV radiatio n intensity at tains 
values in the range E c ~ 3-10 M pc" 2 llSchavdl2004tl . For 
a disc with an exponential surface density profile, he finds: 

R I. j. j.; ix ISC /-. \ 

~faT = 2nh 2 R X c ' ( > 

where M d i sc is the disc mass (gas and stars) and hn is the 
mass scalelength. It is easy to show that when the gas mass 
is negligible this prediction reduces to: 
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Figure 2. K max /fcit versus face-on central surface brightness for 
the V a nd -B-band re sults of Pohlen for the sharply truncated disc 
model iPohlerJl200lT) . The face-on surface brightnesses were cal- 
culated assuming a exponential vertical luminosity distribution. 
As in Fig. lb, the cross-hatched region shows the prediction of 
the collapse model and the dotted line shows the prediction of 
the star-formation treshold model, with the except that here the 
(M/L)+ for the threshold model is arbitrary. 



flmax = ^ (M/L)^Q 

hn E c 

with (M /L)i, the stellar mass-to- light ratio and its central 
surface brightness (linear units). 

The predicted threshold radii according to Eqn. Q arc 
shown in Fig. for three disc masses and Schaye 's (2004) 
fiducial critical surface density E c = 5.9 Mq pc -2 (Ah = 
5.6 10 20 cm -2 ). For the adopted disc masses, the threshold 
model brackets the observations. Interestingly, the model 
predicts an increase in i? max //iR towards small scalelengths. 
Discs with a higher central surface density form stars out to 
a larger radius in terms of scalelengths before reaching the 
critical density. Since, for constant total disc mass, a higher 
surface density disc has a smaller scalelength, it follows that 
smaller scalelength discs have larger Rm^/hn- This anti- 
correlation, which is similar to the obser ved trend in the 
IPohlen et all l|2000h sample ijSchavd^OO^ . is also in accor- 
dance with the present observations. 

Motivated by th e observation that the HI gas fraction in 
the galaxies is small jKreeel. van der Kruit fc de Blokbool 
hereafter KKB), we show Eqn. |2]in Fig. assuming E c = 
5.9 Mqpc -2 and (M/L)* = 2 (constant among galaxies). 
Lower surface brightness discs are predicted to be less ex- 
tended, with a slope in agreement with the observed trend. 
The match is better if we take into account that the inferred 
disc central surface brightnesses are fainter by about 0.7 
magnitudes compared to their face-on counterparts (KKF). 
The scatter and outliers may be explained as being due to 
a non-constancy in E c , e.g. due to a varying UV radiation 
intensity or metallicity, or due to a variation of (M/L)*. 
Altogether, the observations are consistent with a 'critical 
surface brightness' in the range E C /(M/L)* = 1.5-4 Lq { 
pc -2 . Saying it differently, for E c = 5.9 Mqpc -2 the stellar 
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mass-to-light ratio is (M/L)* = 4-1.5 Mq/Lq j. The lower 
part of this range is reasonab le, both from the persp ective of 
stellar population synthesis ijBell fc de Jonell2000F) and ob- 
servations of th e stellar velocity dispersions in galaxy discs 
jBottemalll997h . 

Both models can also be compa red to the truncation 
analysis performed by IPohlenl (1200 ll) in the V and the R 
bands (Fig. |5J . Although obtained using a different method, 
these data show a trend very similar to our /-band result 
(Fig. and to the star-formation treshold model (dotted 
line). Lower surface brightness discs tend to be less extended 
in terms of i?max- 

It has long been known that viscosity driven angular 
momentum re-distribution within a star-forming gas disc 
may dri ve the resulting stel lar disc towards an exponential 
profile <Lin fc Pringlelll987ft . Depending on the details of 
star formation, this process may yield a stellar disc with 
a wel l-defined edge that adva nces radially outward with 
time (iFerguson fc ClarkdlioOlh . Unfortunately, the theory 
has not yet been explicitly investigated with respect to this 
truncation radius. The theory of stoc hastic s elf-propagating 
star formation predicts i? m ax//iR = 4 iSeiden et alJl984F) for 
flat rotation curves. This prediction is similar to that of the 
collapse model, leading to the same difficulties of predicting 
no dependence on surface brightness. 



4 CONCLUSIONS 

At least 20 galaxies in our sample of 34 have truncated stel- 
lar discs, displaying a tight relation between disc scalelength 
and truncation radius. The stellar disc edge seems to occur 
at a larger number of scalelengths in galaxies with a smaller 
scalelength, in agreement with the de crease of -Rmax/^R to - 
wards large scalelengths reported bv IPohlen et alj 1.2000). 
In addition, we observe a tentative correlation between the 
ratio -Rmax/fr-R and the de-projected face-on central surface 
brightness of the disc. These observations appear to have 
two important implications. First, high surface brightness 
spirals with small scalelengths, wh ich are the most numer - 
ous spirals in the local Universe Jde Jong fc Lacevl boOO). 
have an i? ma x//iR of at least four. Secondly, the face-on disc 
surface brightness at the truncation radius is roughly con- 
stant among galaxies. 

The observed truncation radii were compared to the 
predictions of several theories proposed for its origin. The 
data are best reproduced by the c onstant star-formation 
threshold model for the truncation (ISchavd 120041) . In par- 
ticular, this model is able to match the correlation between 
flmax/ftR and surface brightness for reasonable stellar mass- 
to- liglrtj^atioj. TJiecollapse model for disc galaxy formation 
Ivan der Krufdll987l : DSS97), in which the truncation cor- 
responds to the maximum angular momentum in the proto- 
galaxy before disc collapse with detailed conservation of an- 
gular momentum, provides a significantly poorer match to 
the observations. This model would require some amount 
of re-distribution of angular momentum during and/or after 
the collapse. 
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